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CROSS REFERENCE TO RELATED APPLICATIONS 

The present invention claims priority under 35 U.S.C. §119 (e) from U.S. 
provisional application serial no. 60/214781, entitled "Address Generator for LDPC 
10 Encoder and Decoder and Method Thereof," filed June 28, 2000, the contents of which 
are incorporated herein by reference. 

The present invention is a continuation-in-part application of application Serial 
No. 09/559186, entitled "Multi-Mode Iterative Detector" and filed on April 27, 2000, the 
contents of which are incorporated herein by reference, 

15 The present invention is related to the following commonly-assigned, copending 

applications: 

"LDPC Encoder and Method Thereof, filed on even date and assigned 

application Serial No. (Attorney Docket No. MP0064), the contents of 

which are incorporated herein by reference, 

20 "Address Generator for LDPC Encoder and Decoder and Method Thereof 

filed on even date and assigned application Serial No. ___ (Attorney 

Docket No. MP0063), the contents of which are incorporated herein by reference, 
and 

"Parity Check Matrix and Method of Forming Thereof, filed on even date 

25 and assigned application Serial No. (Attorney Docket No. MP0069), 

the contents of which are incorporated herein by reference. 
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BACKGROUND OF THE INVENTION 

Field of the Invention 

The present invention relates generally to a linear block decoder in a data 
transmission system. More particularly, the present invention relates to a low density 
parity-check code (LDPC) decoder a read channel in a disk drive system. 

Description of the Related Art 

Fig. 1 illustrates a conventional digital data transmission system. As shown 
therein, a digital data transmission system comprises a transmitting section 300 for 
transmitting user data to receiver 500 via communication channel 401. 

The operation of transmission section 300 will now be explained. Prior to processing 
by transmitting section 300, input or user data maybe encoded with an error correcting 
code, such as the Reed/Solomon code, or run length limited encoded (RLL) or a 
combination thereof by encoder 302. The encoded output encoder 302 is then 
interleaved by deinterleaver 308 for input to linear block code encoder 304 which 
generates parity data in a known manner utilizing linear block codes. One example of 
a linear block code *is a low-density parity-check code (LDPC) which is discussed by 
Robert G. Gallager in Low-Density Parity-Check Codes, 1963, M.I.T. Press and by 
Zining Wu in Coding and Iterative Detection For Magnetic Recording Channels, 2000, 
Kluwer Academic Publishers, the contents of each of which are incorporated in their 
entirety by reference. Deinterleaver 308 permutes the data so that the same data is 
reordered before encoding, by linear block code encoder 304. By permuting or 
redistributing the data, interleaver 306 attempts to reduce the number of nearest 
neighbors of small distance error events. User data at the output of encoder 302 is 
referred to as being in the channel domain; that is the order in which data is 
transmitted through the channel. The order of data processed by deinterleaver 308 is 
referred to as being in the linear block code domain. The parity data from linear block 
code encoder 304 is combined with the data encoded by encoder 302 by multiplexer 306 
for input to channel transmitter 310. 

Transmitter 310 transmits the combined user and parity data from multiplexer 306 
typically as an analog signal over communication channel 401 in the channel domain. 
Communication channel 401 may include any wireless, wire, optical and the like 
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communication medium. Receiver 500 comprises a front-end circuit 502 comprising 
analog to digital and equalization circuits. The digital signal front-end circuit 502 is 
input to soft channel decoder 504, which provides probability information of the 
detected data. Soft channel decoder may be implemented by a Soft Viterbi Detector or 
5 the like. The output of the soft channel decoder 504, which is in the channel domain, is 
converted into the linear block code domain by deinterleaver 510. Deinterleaver 510 is 
constructed similarly to deinterleaver 308. Soft linear block code decoder 506 utilizes 
this information and the parity bits to decode the received data. One output of soft 
linear block code decoder 506 is fed back to soft channel decoder 504 via interleaver 

10 512, which converts data in the linear block code domain to the channel domain. 
Deinterleaver 512 is constructed to perform the reverse operations of deinterleaver 
510. Soft channel decoder 504 and soft linear block code decoder 506 operate in an 
iterative manner to decode the detected data. 

The other output of soft linear block code decoder 506 is converted from the 

15 linear block domain to the channel domain by interleaver 514. Deinterleaver 514 is 
constructed similarly to interleaver 512. The output of interleaver 514 is passed on for 
further processing to decoder 508. Decoder 508 is implemented to perform the reverse 
operations of encoder 302. 

Summary of the Invention 

20 According to a first aspect of the present invention, a method of decoding low- 

density parity-check codes is provided comprising the steps of (a) calculating llrRmi, for 
each parity check equation, at iteration i-1, in response to step c, (b) decision aided 
equalizing, at iteration i 3 in response to step a; and (c) calculating llrQim, for each 
parity check equation, at iteration i in response to step b. llrQim represents 

25 information from bit node 1 to equation node m, one for each connection, and llrRmi 
represents information from equation node m to bit node 1, one for each connection. 

According to a second aspect of the present invention, a method of decoding low- 
density parity-check codes is provided further comprising the steps of (d) updating, for 
each parity check equation, the smallest llrQim calculated in step c, at iteration i, (e) 
30 updating, for each parity check equation, the second smallest llrQim calculated in step 
c, at iteration i\ and (f) updating, for each parity check equation, an overall sign of 
llrQim, calculated in step c, at iteration L Step a is responsive to step d, step e and 
step f, and step c is responsive to step a. 
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According to a third aspect of the present invention, step a is calculated as 



follows: 



ml 



-s 



i- 1 



sign[llrQ l 1 /,„) • min 1 ' ~ 1 m , if / * i>~ 1 

* mini* otherwise. 



Step c is calculated as follows: 



According to a fourth aspect of the present invention, step c is initialized with 
soft channel information. 

According to a fifth aspect of the present invention, the method further 
comprises the steps of (g) converting | LlrAPPk) information into hard information b c ,k, 
(h) calculating an equation syndrome for each parity check equation, Si. (i) summing 
each of equation syndrome in step h, (j) outputting the hard information if the sum in 
step i is equal to zero; (k) repeating steps a-c if the sum in step i is not equal to zero, (1) 
determining if i is less than a first predetermined value; (m) repeating steps a-c, if in 
step 1, i is less than the first predetermined value; (n) determining if the sum in step i 
is less than a second predetermined value; and (o) outputting the hard information if i 
is at least the first predetermined and the sum in step i is less than the second 
predetermined value. 

According to a sixth aspect of the present invention, the method further 
comprises the steps of (p) for each data k, determining if | llrAPPk | is less than a third 
predetermined value; (q) for each data k, outputting hard information b Ct k if | llrAPPk | 
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is at least the third predetermined value; and (r) for each data k y outputting soft 
channel information b s ,kif | LlrAPPk | is less than the third predetermined value. 

According to a seventh aspect of the present invention, the method further 
comprises the steps of (s) for each data k, determining if a corresponding parity check 
5 equation is violated, (t) for each data k, outputting hard information b c ,k, if the 
corresponding parity check equation is not violated, (u) for each data k, determining if 
| LlrAPPk | is less than a third predetermined value, (v) for each data k, outputting hard 
information b Cl k if | LlrAPPk | is at least the third predetermined value, and (w) for each 
data k, outputting soft channel information b s ,k if | LlrAPPk | is less than the third 
10 predetermined value. 

According to an eighth aspect of the present invention, a computer program for 
decoding low-density parity-check codes is provided comprising the steps of (a) 
calculating llrRmi, for each parity check equation, at iteration in response to step c, 
(b) decision aided equalizing, at iteration i, in response to step a; and (c) calculating 
15 llrQim, for each parity check equation, at iteration i in response to step b. llrQim 
represents information from bit node 1 to equation node m, one for each connection, 
and LlrRmi represents information from equation node m to bit node 1, one for each 
connection. 

According to a ninth aspect of the present invention, a computer program for 
20 decoding low-density parity-check codes is provided further comprising the steps* of (d) 
updating, for each parity, check equation, the smallest llrQim calculated in step c, at 
iteration i } (e) updating, for each parity check equation, the second smallest llrQim 
calculated in step c, at iteration and (f) updating, for each parity check equation, an 
overall sign of llrQim, calculated in step c, at iteration i. Step a is responsive to step d, 
25 step e and step f, and step c is responsive to step a. 

According to a tenth aspect of the present invention, step a is calculated as 
follows: 
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Step c is calculated as follows: 



According to an eleventh aspect of the present invention, step c is initialized 
with soft channel information. 

According to a twelfth aspect of the present invention, the computer program 
further comprises the steps of '(g) converting | LLrAPPk | information into hard 
information b c ,k, (h) calculating an equation syndrome for each parity check equation, 
Si. (i) summing each of equation syndrome in step h, (j) outputting the hard 
information if the sum in step i is equal to zero; (k) repeating steps a-c if the sum in 
step i is not equal to zero, (1) determining if i is less than a first predetermined value; 
(m) repeating steps a-c, if in step 1, i is less than the first predetermined value; (n) 
determining if the sum in step i is less than a second predetermined value; and (o) 
outputting the hard information if i is at least the first predetermined and the sum in 
step i is less than the second predetermined value. 

According to a thirteenth aspect of the present invention, the computer program 
further comprises the steps of (p) for each data k, determining if | LLrAPPk | is less than 
a third predetermined value; (q) for each data k, outputting hard information b c ,k if 
| LLrAPPk | is at least the third predetermined value; and (r) for each data k, outputting 
soft channel information b s ^ if | LLrAPPk | is less than the third predetermined value. 

According to a fourteenth aspect of the present invention, the computer program 
further comprises the steps of (s) for each data k, determining if a corresponding parity 
check equation is violated, (t) for each data k, outputting hard information b c ,k, if the 
corresponding parity check equation is not violated, (u) for each data k, determining if 
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| LlrAPPk | is less than a third predetermined value, (v) for each data k, outputting hard 
information b c ,k if | LlrAPPk | is at least the third predetermined value, and (w) for each 
data k, outputting soft channel information b s ,k if | LlrAPPk | is less than the third 
predetermined value. 

According to a fifteenth aspect of the present invention, a decoder for decoding 
low-density parity-check codes is provided comprising first calculating means for 
calculating LlrRmi, for each parity check equation, at iteration M, decision aided 
equalizing means for equalizing LLrRmi, at iteration i, in response to the first 
calculating means, and second calculating means for calculating LLrQun, for each parity 
check equation, at iteration i in response to the decision aided equalizing means. 
LLrQim represents information from bit node 1 to equation node m, one for each 
connection, and LLrRmi represents information from equation node m to bit node 1, one 
for each connection, and the first calculating means is responsive to the second 
calculating means. 

According to a sixteenth aspect of the present invention, The decoder, further 
comprises memory means for storing for each parity check equation, the smallest 
llrQim calculated by the second calculating means, at iteration i, the second smallest 
llrQim calculated by the second calculating means, at iteration i\ and an overall sign of 
llrQim, calculated by the second calculating means, at iteration L The first calculating 
means is responsive to the memory means, and the second calculating means is 
responsive to the first calculating means. 

According to a seventeenth aspect of the present invention, the first calculating 
means calculates as follows: 



UrR 
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The second calculating means calculates as follows: 
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According to an eighteenth aspect of the present invention, the second 
calculating means is initialized with soft channel information. 

According to a nineteenth aspect of the present, invention, the decoder, further 
5 comprises an address generating means for providing indices of each of the parity 
check equations. 

According to a twentieth aspect of the present invention, the second calculating 
means provides an index of the smallest UrQim. 

According to a twenty-first aspect of the present invention, multiplying means is 
10 provided for scaling UrR m i, calculated by the first calculating means. 

According to a twenty-second aspect of the present invention, the decoder 
further comprising slicing means for converting | LLrAPPk | information into hard 
information b c ,k; equation vector means for calculating an equation syndrome for each 
parity check -equation Si, summing means for summing each equation syndrome 
15 calculated by the equation vector means. The hard information is output if the sum by 
the summing means is equal to zero, and the calculations by the first and second 
calculating means are repeated if the sum summed by the summing means is not equal 
to zero. 

According to a twenty-third aspect of the present invention, the decoder further 
20 comprises first threshold means for determining if i is less than a first predetermined 
value, and second threshold means for determining if the sum by the summing means 
is less than a second predetermined value. The calculations by the first and second 
calculating means are repeated if i is less than the first predetermined value as 
determined by the first threshold means, and the hard information is output if i is at 
25 least the first predetermined and the sum by the summing means is less than the 
second predetermined value. 

According to a twenty-fourth aspect of the present invention, the decoder further 
comprises third threshold means for determining for each data k if | LLrAPPk | is less 
than a third predetermined value. For each data k, hard information b c ,k is output if 
30 | LLrAPPk | is at least the third predetermined value as determined by the third 
threshold means; and for each data k, soft channel information b St k is output if 
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| LLr APPk | is less than the third predetermined value as determined by the third 
threshold means. 

According to a twenty-fifth aspect of the present invention, the decoder further 
comprises judging means for determining for each data k if a corresponding parity 
5 check equation is violated and third threshold means for determining for each data k if 
| llrAPPk | is less than a third predetermined value. For each data k, hard information 
b c ,k, is output if the corresponding parity check equation is not violated as determined 
by the judging means, for each data k, hard information b c ,k is output if | llrAPPk | is at 
least the third predetermined value as determined by the third threshold means; and 
10 wherein for each data k, soft channel information b s jc is output if | llrAPPk j is less than 
the third predetermined value as determined by the third threshold means. 

According to a twenty-sixth aspect of the present invention, a decoder for 
decoding low-density parity-check codes is provided comprising a first calculator to 
calculate UrRmi, for each parity check equation, at iteration a decision aided 

15 equalizer to equalize UrRmi, at iteration i, in response to the first calculator; and a 
second calculator to calculate UrQim, for each parity check equation, at iteration i in 
response to the decision aided equalizer. UrQim represents information from bit node 1 
to equation node m, one for each connection, and UrRmi represents information from 
equation node m to bit node 1, one for each connection. The first calculator is 

20 responsive to the second calculator. 

According to a twenty-seventh aspect of the present invention, the decoder 
further comprises a memory to store for each parity check equation, the smallest 
UrQim calculated by the second calculator, at iteration i; the second smallest UrQim 
calculated by the second calculator, at iteration i; and an overall sign of UrQim, 
25 calculated by the second calculator, at iteration L The first calculator is responsive to 
the memory, and the second calculator is responsive to the first calculator. 

According to a twenty-eighth aspect of the present invention, the first calculator 
calculates as follows: 
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The second calculator calculates as follows: 



nrQj- UrP\+ £ UrRj' 1 

ir;' * *rt 

According to a twenty-ninth aspect of the present invention, the second 
5 • calculator is initialized with soft channel information. 

According to a thirtieth aspect of the present invention, the decoder further 
comprises an address generator for providing indices of each of the parity check 
equations. 

According to a thirty-first aspect of the present invention, the second calculator 
10 provides an index of the smallest ItrQim. 

According to a thirty-second aspect of the present invention, the decoder further 
comprises a multiplier to scale LtrRmi, calculated by the first calculator. 

According to a thirty-third aspect of the present invention, the decoder further 
comprises a slicer to convert | LlrAPPk | information into hard information b c ,k ;j and an 
15 equation vector circuit to calculate an equation syndrome for each parity check 
equation, Si. 

According to a thirty-fourth aspect of the present invention, the decoder further 
comprises a summer to sum each equation syndrome calculated by the equation vector 
circuit. The hard information is output if the sum summed by the summer i is equal to 
20 zero, and the calculations by the first and second calculator are repeated if the sum 
summed by the summer is not equal to zero. 

According to a thirty-fifth aspect of the present invention, the decoder further 
comprises a first threshold detector to determine if i is less than a first predetermined 
value; and a second threshold detector to determine if the sum by the summer is less 
25 than a second predetermined value. The calculations by the first and second calculator 



UrR 1 , 1 = < 
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are repeated if i is less than the first predetermined value as determined by the first 
threshold detector; and the hard information is output if i is at least the first 
predetermined and the sum by the summing means is less than the second 
predetermined value. 

5 According to a thirty-sixth aspect of the present invention, the decoder further 

comprises a third threshold detector to determine for each data k if | LlrAPPk | is less 
than a third predetermined value. For each data k, hard information b c ,k is output if 
| LLrAPPk | is at least the third predetermined value as determined by the third 
threshold detector; and for each data k } soft channel information b Sj k is output if 
10 | LLrAPPk | is less than the third predetermined value as determined by the third 
threshold detector. 

According to a thirty-seventh aspect of the present invention, the decoder 
further comprises a judging circuit to determine for each data k if a corresponding 
parity check equation is violated, and a third threshold detector to determine for each 

15 data k if | LLrAPPk | is less than a third predetermined value. For each data k, hard 
information be* is output, if the corresponding parity check equation is not violated as 
determined by the judging circuit, and for each data k, outputting hard information b c ,k 
if | LLrAPPk | is at least the third predetermined value as determined by the third 
threshold detector; and for each data k, soft channel information b s ,k is output if 

20 | LLrAPPk | is less than the third predetermined value as determined by the third 
threshold detector. 

Other objects and attainments together with a fuller understanding of the 
invention will become apparent and appreciated by referring to the following 
description and claims taken in conjunction with the accompanying drawings. 

25 Brief Description of the Drawings 

In the drawings wherein like reference symbols refer to like parts. 

Fig. 1 is a block diagram of a conventional data transmission system; 

Fig. 2 is a block diagram of a data transmission system in accordance with the 
present invention; 

30 Fig. 3 is . an example of a parity check matrix in accordance with the present 

invention; 
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Fig. 4 is a block diagram of a memory arrangement of the decoder in accordance 
with the present invention; 

Fig. 5 is a block diagram of the decoder in accordance with the present 
invention; 

Fig. 6 is an example of a factor graph; 

Fig. 7 is an example of a factor graph illustrating the calculation of LlrQim; 

Fig. 8 is an example of a factor graph illustrating the calculation of LLrRmi; 

Fig. 9 is a graph of a function for calculating a soft decision; 

Fig. 10 is a block diagram of a decision circuit of the decoder of Fig. 5; 

Fig. 11 is a block diagram of an alternate decision circuit of the decoder of Fig. 5; 

Fig. 12 is a factor graph of the decoding procedure in accordance with the 
present invention; 

Fig. 13 is a flow diagram of the decoding procedure in accordance with the 
present invention; 

Fig. 14 is a flow chart of the decision procedure of the decision circuit of Fig. 10; 

Fig. 15 is a flow chart of the decision procedure of the alternate decision circuit 
of Fig. 11; 

Fig. 16 is a block diagram of a conventional decoding implementation; and 

Fig. 17 is a block diagram of a decoding implementation in accordance with the 
present invention. 

Description of the Preferred Embodiments 

Reference is now made to Fig. 2, which is a block diagram of a data transmission 
system in accordance with the present invention. In general as shown therein, a 
digital data transmission system comprises a transmitting section 300' for 
transmitting user data to receiver 500' via communication channel 401- The inventors 
have observed that a linear block code encoder is not dependent on a position of a bit 
interleaved. Rather the linear block code encoder only requires a list of equations for a 
given bit. In other words, there is no need to process the data in the order defined by 
the interleave^, instead data may be processed in the same order as it is written to the 
channel. This can be accomplished by incorporating an address generator to provide 
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an address of the appropriate equation of the linear block code encoder. This principle 
can be similarly applied to the soft linear block decoder. As a result, deinterleaver 308 
of the conventional system is now replaced by address generator 328, and 
deinterleaver 510 is now replaced by address generator 530. Accordingly, there is no 
5 requirement for the physical interleaving of data in the receiver 500', since the data 
remains in the same order as the order of bits of data in the channel throughout this 
system. The order of bits of data transmitted through the channel is referred to as the 
channel domain. 

The operation of transmission section 300* will now be explained. Prior to 
10 processing by transmitting section 300', as in the conventional system, input or user 
data maybe encoded with an error correcting code, such as the Reed/Solomon code, or 
run length limited encoded (RLL) or a combination thereof by encoder 302. Addresses 
for the parity equations of linear block code encoder 304 are generated by address 
generator 328 in accordance with an index of the bits of data, the index being 
15 determined by address generator 328. Address generator 328 is responsive to counter 
730 under the control of controller 740. Controller 740 synchronizes counter 730 to the 
output of encoder 302 so that counter 730 can provide a count of the number of bits in a 
codeword output by encoder 302 and a count of the number of codewords. In the 
preferred embodiment the data block size is 5000 bits. 

20 Referring again to Fig. 2, linear block code encoder .304 utilizes the user data 
and address from address generator 328 to provide the parity bits to multiplexer 306. 
Address generator 328 is described in more detail in commonly assigned, copending 
patent application entitled "Address Generator for LDPC Encoder and Decoder and 
Method Thereof 5 filed on even date and assigned application Serial No. 

25 (Attorney Docket No. MP0063), the contents of which are incorporated herein by 
reference. Linear block code encoder 304 is preferably implemented as a low-density 
parity-check code (LDPC) encoder as described in commonly assigned, copending 
patent application entitled "LDPC Encoder and Method Thereof, filed on even date 
and assigned application Serial No. (Attorney Docket No. MP0064), the 

30 contents of which are incorporated herein by reference. The parity data from linear 
block code encoder 304 is combined with the data encoded by encoder 302 by 
multiplexer 306 for input to channel transmitter 310. In the preferred embodiment, the 
combined data consists of series of a pair parity bits followed by 34 bits of user data. 
This constraint is established by RLL constraint encoder 302. 
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Transmitter 310 transmits the combined user and parity data from multiplexer 
306 typically as an analog signal over communication channel 401 in the channel 
domain: Communication channel 401 may include any wireless, wire, optical, 
magnetic and the like. 

Receiver 500' comprises an analog to front-end circuit 502 comprising analog to 
digital and equalization circuits. The digital signal from front-end circuit 502 is input 
to soft channel decoder 504, which provides soft or probabilistic information of the 
detected data to soft linear block decoder 506. Soft channel decoder may be 
implemented as a Soft Viterbi Detector or the like, and address generator 530 may be 
constructed similarly as address generator 328 in transmission section 300'. In the 
preferred the Soft Viterbi Detector may be implemented with a Soft-Output Viterbi 
Algorithm which is described in J. Hagenauer and P. Hoeher: "A Viterbi algorithm 
with soft-decision outputs and its applications", Proc. IEEE GLOBECOM '90, Dallas, 
Tex., pp. 47.1.1-47.1.7, November 1989, the contents of which are incorporated by 
reference. 

The soft information output by soft channel decoder 504 remains in the channel 
domain and is decoded by soft linear block code decoder 506, in accordance with the 
address of the parity equations generated by address generator 530. Address generator 
530 is responsive to counter 735 under the control of controller 745. Controller 745 
synchronizes counter 735 to the output of soft channel decoder 504 so that counter 735 
can provide a count of the number of bits in a codeword output by soft channel decoder 
504 and a count of the number of codewords. 

Soft linear block code decoder 506 operates in combination with soft channel 
decoder 504 and address generator 530 in an iterative fashion. Soft linear block code 
decoder is preferably implemented as a low-density parity- check code (LDPC) decoder, 
as described in detail hereinbelow. It is noted that since the soft information from soft 
channel decoder 504 to soft linear block code decoder 506 are both in the channel 
domain, thus as noted above, there is no need for any interleaves or deinterle avers in 
receiver 500'. 

After the iterative process has completed, the output of soft linear block code 
decoder 506 is passed on for further processing to decoder 508. Decoder 508 is 
implemented to perform the reverse operations of encoder 302 or correct for any data 
errors. 
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Prior to discussing the construction and operation of the LPDC decoder, 
reference is now made to Fig. 3 for an explanation of the parity check matrix. Fig. 3 
shows only one example of a parity check matrix. Of course, other parity check 
matrices are contemplated. The preferred matrix is 222 rows (or equations) by 5402 
5 columns, which comprises 220 linearly independent rows (where 5402 =73*74). The 
matrix can be divided into three tiers of equations having 73, 74 and 75 equations, 
respectively. The set of independent rows can be obtained by canceling the last row of 
the second tier and third tier, namely the 147 th row and the 222 nd row. As shown in 
Fig. 3, the following table shows the values of the elements in the matrix: 



Tier 


1 th position 


I th position 


1 


lifr = i(mod73) 


0ifr*i(mod73) 


2 


lifr = i(mod74) 


0 if r* i(mod74) 


3 


lifr = i(mod75) 


0ifr*i(mod75) 



10 

A matrix having 5402 columns can process a maximum LDPC codeword of 5402 
bits. Of course,, as will be appreciated by one of ordinary skill in the art, the matrix 
may be truncated to accommodate a smaller block, however the matrix must be at least 
222 x 4366 which is dependent on the constraint of encoder 302. This constraint is for 

15 example a RLL constraint. The preferred matrix contains no cycles, since a matrix 
having cycles has degraded performance that degrades significantly. With the first 
tier only, the parity check matrix has a Dmin=2; by adding the second tier, the parity 
check matrix has a Dmin=4; and by adding the third tier, the parity check matrix has a 
Dmin=6. A further description of the parity check matrix is provided in commonly 

20 assigned, copending application entitled, "Parity Check Matrix and Method of 

Designing Thereof, filed on even date and assigned application Serial No. 

(Attorney Docket No. MP0069), the contents of which are incorporated herein by 
reference. 

The LDPC decoder is preferably implemented using the sum-product algorithm, 
25 which is described by Zining Wu in Coding and Iterative Detection For Magnetic 
Recording Channels, 2000, Kluwer Academic Publishers, the contents of each of which 
are incorporated in its entirety by reference. 

A linear code is a set of codewords, x, satisfying the matrix equation (1) 
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H«x=0 (1), 
where H is an M x L matrix, and x is a 1 x L vector. 

The parity check matrix for an LDPC encoder/decoder is sparse, that is a small 
portion of the elements being one, all other elements being zero. An example of a 
parity check matrix is shown in equation 2 that corresponds to the factor graph shown 
in Fig. 6. 



H = 



1011 
0111 



(2) 



which defines the following parity-check equations for the codeword x 



xl+x3+x4=0 
x2+x3+x4=0 

Fig. 6 is a factor graph of the parity matrix of equation 2. The factor graph 
contains two types of nodes the bit node (e.g.. bl, b2, b3, and b4) and the check nodes 
(e.g. el, e2). Each bit node corresponds to a bit in the codeword, and each check node 
represents a parity-check equation (i.e., a row in the parity check matrix H). Hence, 
the factor graph for an LDPC code with an M x L parity check matrix H contains M 
check nodes and L bit nodes. An edge between a check node and a bit node exists if 
and only if the bit participates in the parity-check equation represented by the check 
node. The factor graph shown in Fig. 3 is "bipartite" in which the nodes can be 
separated to two groups, namely check nodes and bit nodes. Connections are allowed 
only between nodes in different groups. 

A cycle in a factor graph refers to a finite set of connected edges that start and 
end at the same node. The bold lines in Fig. 6 represent a cycle length of four. As can 
be appreciated by one of ordinary skill in the art, four is the shortest cycle length a 
parity check matrix can have. 

In the sum-product decoding algorithm, there are four types of soft information 
circulating in the decoding process: 

1. ILrPi: The soft information from channel APP (a posteriori probability) 
decoder or soft channel decoder 504 for the first iteration or from decision aided 
equalization circuit 856 for subsequent iterations, one for each bit. 
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2. LlrQun: The information from bit node 1 to equation node m, one for each 
connection. 

3. llrRmi: The information from equation node m to bit node 1, one for each 
connection. 

5 4. llrAPPi: the overall soft information after each iteration, one for each bit. 

Each iteration performs three steps of operations (superscripts denote iteration 
number): 

1. Each bit calculates the information 11 rQ that it passes to the 
connecting equations, which is the sum of the extrinsic information llrR from 
10 the previous iteration with the channel information tlrP from the current 

iteration, excluding the extrinsic information from the same equation, i.e.: 

LlrQim' = LLtPl' + X LlrRm'i'" 1 (3) 

2 Each check equation (each row) calculates the extrinsic information 

15 for each involved bit. The "LLRXOR" (2 ©) denotes the LLR operation discussed 

below. 



llrfW = Z © llrQi'm 1 (4) . 

r*i 

20 3. Each bit calculates the overall soft information by summing up all 

the UrRmi's and the llrPi 

llrAPPi': = UrPi' +2 llrfW (5) 
m 

The "LLRXOR" operation is defined as: 
25 LLRXOR(y,z) = log(ey+e z )-log(ey«+e°) (6) 

max(y,z)+log(l+e- 1 y- 2 i )-max(x+y,0)-log(l+e- ly +z I ) 

where log (l+e-l°l) can be implemented with a look-up table. Such an operation 
. requires 5 additions, 2 absolute values and two table look-ups from a look-up table. 
30 This procedure is referred to as the "Sum Product" (SP) algorithm. 
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The "LLRXOR" operation can be further simplified by eliminating the table look- 
ups, i.e. omitting the log (l+e-l*!) term in equation (6). The equation is simplified as 
follows: 

LLRXOR(y,z)=-sign(y)sign(z)min( | y | , | z | ) (7). 
5 As is apparent to one of ordinary skill in the art, equation (7) requires only 

absolute value and comparison operations. This operation is referred to as the "Signed 
Sum Product" (SSP) algorithm. 

In the SSP algorithm, the magnitudes of the extrinsic information conveyed in 
each equation are determined by the two log-likelihood ratios with the smallest 
10 amplitudes. Other LLR's only contribute to the sign of the information, not to the 
amplitude. Fig. 7 shows the LlrCWs passed from bits to an equation, and Fig. 8 shows, 
the LlrRmi's relayed by the equation. 

A conventional implementation of SP or SSP requires L units of memory to store 
the UrPi:'s; L units of memory to store the UrAPPi/s; ML units of memory to store the 

15 . LlrCWs; and ML units of memory to store the LLrfWs. Soft information is updated bit 
by bit within each iteration. These implementations are referred to as the "bit-based" 
implementations. As will be appreciated by one of ordinary skill in the art, the bit- 
based implementations require a large amount of memory and increased hardware 
complexity. Fig. 12 illustrates an implementation of SSP, which is generally 

20 performed in a serial manner. As shown therein, soft channel decoding is performed 
for the first iteration. Thereafter ILrQim and LLrRmi are sequentially calculated. For 
subsequent iterations i, decision aided equalization is then performed on the 
information. After the decision aided equalization iteration LlrQim is calculated in 
accordance with equation (3), and then LLrRmi- is calculated in accordance with 

25 equation (4). The conventional process requires, for a 4000 bit codeword, 
approximately 4000 cycles for soft channel decoding, approximately 4000 cycles for 
calculating LlrQim, approximately 4000 cycles for calculating LLrRmi, and approximately 
4000 cycles for performing decision aided equalization. The amount of time required 
by conventional methods degrades performance. 

30 In accordance with the present invention, a pipelined architecture is proposed 

for the SSP circuit, as shown in Fig. 17. In accordance with the present invention, soft 
channel decoding is performed to initialize the SSP circuit, and for iteration i, LLrQim 1 is 
calculated. At the same time LLrRmi 1 " 1 for iteration i-1 is also calculated. Based on this 
result and LLrAPPi decision aided equalization is then performed. As can be 

35 appreciated by one of ordinary skill in the art, the pipelined architecture allows for 
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reduced pendency, resulting in higher throughput In accordance with the present 
invention, each iteration takes about 4000 cycles for a 4000 bit codeword. According to 
another aspect of the present invention, an "equation-based" arrangement is proposed 
which implements the same algorithm with much less hardware complexity as 
5 described herein below. 

The "equation-based" arrangement takes advantage of both the sparseness of 
the parity-check matrix and the simplicity in the SSP decoding. It is noted that the 
number of equations for each bit involved corresponds to the column weight of the 
parity check matrix, which in the preferred embodiment is three. Of course other 

10 appropriate column weights may be employed. According to equation (7), each 
equation node only passes two values to the connected bit nodes if signs of the soft 
information are not considered. Therefore, there is only a need to store two LLR 
magnitudes for each equation, and the signs for all the UrQtm's. for each equation, the 
two smallest LLR's (mini and min2), the position of mini (i), and the sign of the 

15 equation(s); for each bit are stored. Additionally, three signs are stored corresponding 
to the signs of three UrQim's from the previously iteration. UrRim's can be calculated 
from the stored information as follows: 



UrR* , ' ' = < 
ml 



-s* f ~ * ■ signal rO* '/mj * mini 1 ~ otherwise. 



20 where SmM-1) IT -signdlrCW) is the sign of the mth equation 

1 

(8) ' 

Referring to block diagram in Fig. 5, the corresponding flow chart in Fig. 13, 
and the corresponding factor graph in Fig. 12, the "equation-based" SSP circuit 
25 performs the following process: 

1. The equation-based SSP circuit is initialized with the output of soft 
channel decoder 504 for the calculation of LLrPt, which is selected by multiplexer 854. 
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During subsequent iterations, multiplexer 854 selects the output of decision aided 
equalization circuit 856, as discussed hereinbelow 

2. At each time clock, address generator 804 calculates the indices of each 
equation for the current bit, and position index generator 802 calculates the position 

5 index thereof. Three llrfWs are calculated using Equation 8 by a calculating means or 
calculator or equation 1 LLR update circuit 816, equation 2 LLR update circuit 818, 
and equation 3 LLR update circuit 820. 

3. Each bit calculates the information that it passes to the connecting 
equations by calculating means or calculator as follows, 



circuits 832, 834, and 842, respectively, to synchronize with llrPi. 

4. Each bit stores the three signs of the LlrQ's in memory means or memory 
15 sections 708, 712 and 726 for equations 1, 2 and 3 respectively. 

5. Each bit updates the equation information (mini, min2, i, and s) for the 
three equations involving the current bit by equation 1 update circuit 808, equation 2 
update circuit 810, and equation 3 update circuit 830, respectively, in accordance with 
the following procedure: 



10 



(9) 



The UrRmi's are multiplied by a scaling factor by multipliers 836, 838 and 840 
respectively to reduce circulation of any errors. The UrRmi's are also delayed by delay 



20 




minV m ^minV m 



else if 



25 



endif 
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5. Each bit calculates the overall soft information by summing up all the 
LLrRmi's, and llrPi (via delay circuit 848) as an input to decision aided equalizer 856. 



m 

5 Log-likelihood ratios are converted to probability information by converter 850. 

Converter 850 implements the function shown in Fig. 9. In the preferred embodiment, 
the function is implemented as a piecewise linear approximation, shown as the dashed 
lines in Fig. 9. Decision aided equalizer 856 processes llrAPPi and samples stored in 
samples memory 864 from an equalizer (not shown) in front end 502. Decision aided 
10 equalization is described in more detail by Zining Wu in Coding and Iterative Detection 
For Magnetic Recording Channels, 2000, Kluwer Academic Publishers, the contents of 
which are incorporated in its entirety by reference. 

Decision circuit 900, which is responsive to the output of summer 852, address 
generator 530 and soft channel decoder 504 determines when to terminate this 
15 iterative process, as is described in detail below. 

Fig. 4 is a block diagram of a memory arrangement of the decoder in accordance 
with the present invention. As shown therein, memory or memory means 750 has a 
partition for each equation, and each equation has a partition for storing magnitude 
data, index data and sign data for iteration i and for iteration As used herein, the 

20 term partition refers either to a . physical memory partition or a logical memory 
partition. Equation 1 update circuit 808, for each clock cycle, writes during iteration i 
to memory partition 806a or memory partition 806b i, and reads data previously 
stored from iteration i-1 from memory partition 806b or memory partition 806a. 
Multiplexer 902 selects the memory partition for write operations and multiplexer 904 

25 selects the memory partition for read operations. Equation 1 LLR update circuit 816 
reads data previously stored from iteration i-1 from memory partition 806b or memory 
partition 806a as selected by multiplexer 914. Equation 2 update circuit 810, equation 
3 update circuit 830, equation 2 LLR update circuit 818, and equation 3 LLR update 
circuit 820 operate similarly. 

30 Fig. 10 is a block diagram of decision circuit 900, and Fig. 14 is a flow chart 

thereof. Decision circuit 900 is responsive to the output from soft channel decoder 914, 
llrAPPi from summer or summing means 852 and address generator or address 
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generator means 530. Decision circuit 900 comprises memory 914 for storing hard 
decision data b s ,k (for k= 1:5402) from soft channel decoder 504, sheer or slicing means 
888 for converting | tlrAPPi | information into hard information b c ,k, memory 904 for 
storing b C) k, a threshold detector 902 to determine if | LlrAPPi | from summer 852 is less 
5 then a threshold value (t2). Decision circuit or decision means 900 further comprises 
counter 920 for counting the number of iterations performed by the decoder and 
threshold detector 910 to determine if the number of iterations exceeds threshold value 
t3, which in the preferred embodiment is four. Equation vector circuit equation vector 
means 908 determines the equation syndromes, one for each equation, Sm, m=l:222 in 
10' which S m =0 when equation i is satisfied or S m =l when equation m is violated. 
Summer 906 sums each of the Sm m=l:222. 

The operation of decision circuit 900 will now be explained in conjunction with 
the flow chart of Fig. 14. After an iteration of iterative decoding, the equation 
syndrome is calculated by equal vector circuit 908 and summed by summer or 

15 summing means 906. Threshold circuit 912 determines if the sum of all S m =0, then a 
signal is sent to multiplexer 854 to output b c ,k from memory 904. If the sum of all S m 
*0, then a determination is made by threshold 910 whether the number of iterations is 
less than threshold value t3. If the number of iterations is less than threshold value t3, 
the next iteration is performed. Otherwise threshold circuit 912 determines if the sum 

20 of all S m is less than threshold ti. If so, then a signal is sent to multiplexer 854 to 
output b C) k from memory 904. Alternatively, the next step in the process is to 
determine if |UrAPPi| is less then threshold value (t2>, if so, multiplexer 854 selects 
b s ,k to be output. If | LlrAPPi | is not less then threshold value (t2), multiplexer 854 
selects b c ,k to be output. 

25 Fig. 11 is a block diagram of decision circuit 900' and Fig. 15 is a flow chart 

thereof. Decision circuit 900' is responsive to the output from soft channel decoder 914, 
LLrAPPi from summer 852 and address generator 530. Decision circuit 900' comprises 
memory 914 for storing hard decision data b s ,k (for k=l:5402) from soft channel decoder 
504, slicer 888 for converting | LLrAPPi | information into hard information b c ,k, memory 

30 904 for storing b Ct k, a threshold detector 902 to determine if | LLrAPPi | from summer 
852 is less then a threshold value (t2). Decision circuit 900' further comprises counter 
920 for counting the number of iterations performed by the decoder and threshold 
detector 910 to determine if the number of iterations exceeds threshold value t3, which 
in . the preferred embodiment is four. Equation vector circuit 908 determines the 

35 equation syndromes, one for each equation, S m , i=l:222 in which S m =0 when equation i 
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is satisfied or S m =1 when equation i is violated. Summer 906 sums each of the S m 
m=l:222. 

The operation of decision circuit 900' will now be explained in conjunction with 
the flow chart of Fig. 15. After an iteration of iterative decoding, the equation 
5 syndrome is calculated and summed by equal vector circuit 908 and summer 906. 
• Threshold circuit 912 determines if the sum of all S m =0, then a signal is sent to 
multiplexer 854 to output b c ,k from memory 904. If the sum of all S m *0, then a 
determination is made by threshold 910 whether the number of iterations is less than 
threshold value t3. If the number of iterations is less than threshold value t3, the next 

10 iteration is performed. Otherwise threshold circuit 912 determines if the sum of all S m 
is less than threshold ti. If so, then a signal is sent to multiplexer 854 to output b c ,k 
from memory 904. Alternatively, the next step in the process is to determine if bk is 
involved in an equation that is violated. If bk is not involved in an equation that is 
violated multiplexer 854 b c ,k is provided as the output. Otherwise threshold 902 

15 determines if | llrAPPk | is less then threshold value Oa), if so multiplexer 854 selects 
b s ,k to be output. If | llrAPPk | is not less then threshold value (t2), multiplexer 854 
selects b C) k to be output. 

While the invention has been described in conjunction with several specific 
embodiments, it is evident to those skilled in the art that many further alternatives, 

20 modifications and variations will be apparent in light of the foregoing description. 
More specifically, while the present - invention is preferably implemented as an 
integrated circuit, it is contemplated that the present invention may also be 
implemented as discrete components or a general-purpose processor operated in 
accordance with program code instructions or computer program or combination 

25 thereof. These program code instructions can be obtain from a medium, such as 
network, local area network, the Internet, or storage devices. Such storage devices 
include, by way of example, magnetic storage devices, optical storage devices, 
electronic storage devices, magneto-optical device, and the like. Thus, the invention 
described herein is intended to embrace all such alternatives, modifications, 

30 applications and variations as may fall within the spirit and scope of the appended 
claims. 
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WHAT IS CLAIMED IS: 

1. A method of decoding low-density parity -check codes comprising the steps of: 

a. calculating llrRmi, for each parity check equation, at iteration in 
response to step c; 

b. decision aided equalizing, at iteration i, in response to step a; and 

c. calculating llrQim, for each parity check equation, at iteration i in 
response to step b, 

wherein llrQim represents information from bit node I to equation node m, one 
for each connection, and 

wherein llrRmi represents information from equation node m to bit node 1, one 
for each connection. 

2. The method of Claim 1, further comprising the steps of: 

d. updating, for each parity check equation, the smallest llrQim calculated in 
step c, at iteration z; 

e. updating, for each parity check equation, the second smallest llrQim 
calculated in step c, at iteration i\ and 

f. updating, for each parity check equation, an overall sign of llrQim, 
calculated in step c, at iteration i, 

wherein step a is responsive to step d, step e and step f, and 

wherein step c is responsive to step a. 

3. The method of Claim 2, wherein step a is calculated as follows: 
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ml 



-s i ~ ] 
m 



sign(llrO' 1 /„,) min l' ~ 1 m /\fl * l'^ 1 
^~ 1 • signer O' ~ '/^j m//i2 /_ othe 



rwise. 



4. The method of Claim 2, wherein step c is calculated as follows: 



HrQj-llrP}* T UrRj- 1 

to' * «t 



5. The method of Claim 2, wherein step c is initialized with soft channel 
5 information. 

6. The method of Claim 2, further comprising the step of providing indices of each 
of the parity check equations. 

7. The method of Claim 2, further comprising the step of providing an index of the 
smallest llrQim. 

10 8. The method of Claim 2, further comprising the step of scaling UrR m i, calculated 
in step a. 

9. The method of Claim 52, further comprising the steps of: 
i summing each of equation syndrome in step h, 

j outputting the hard information if the sum in step i is equal to zero; and 
15 k repeating steps a-c if the sum in step i is not equal to zero. 
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10. The method of Claim 9 3 further comprising the steps of: 

1 determining if i is less than a first predetermined value; 

m repeating steps a-c, if in step 1, i is less than the first predetermined 

value; 

5 n determining if the sum in step i is less than a second predetermined 

value; and 

o outputting the hard information if i is at least the first predetermined and 
the sum in step i is less than the second predetermined value. 

11. The method of Claim 10, further comprising the steps of: 

10 p for each data k, determining if | LLrAPPk | is less than a third 

predetermined value; 

q for each data k, outputting hard information b c ,k if | LLrAPPk | is at least 
the third predetermined value; and 

r for each data k, outputting soft channel information b s ,k if | LLrAPPk | is 
15 less than the third predetermined value. 

12. The method of Claim 10, further comprising the steps of: 

s for each data k, determining if a corresponding parity check equation is 
violated; 

t for each data k, outputting hard information b c ,k, if the corresponding 
20 parity check equation is not violated; 

u for each data k } determining if | LLrAPPk | is less than a third 
predetermined value; 

v for each data k, outputting hard information b c x if | LLrAPPk | is at least 
the third predetermined value; and 
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w for each data k, outputting soft channel information b s ,k if | LLrAPPk | is 
less than the third predetermined value. 

13. A computer program for decoding low-density parity-check codes comprising the 
steps of: 

5 g. calculating llrRmi, for each parity check equation, at iteration M, in 

response to step c; 

h. decision aided equalizing, at iteration i, in response to step a; and 

i. calculating llrQim, for each parity check equation, at iteration i in 
response to step b; 

10 wherein llrQim represents information from bit node 1 to equation node m, one 

for each connection, and 

wherein UrR m i represents information from equation node m to bit node 1, one 
for each connection. 

14. The computer program of Claim 13, further comprising the steps of: 

15 j. updating, for each parity check equation, the smallest llrQim calculated in 

step c,.at iteration i\ 

k. updating, for each parity check equation, the second smallest llrQim 
calculated in step c, at iteration i; and 

1. updating, for each parity check equation, an overall sign of llrQim, 
20 calculated in step c, at iteration /, 

wherein step a is responsive to step d, step e and step f, and 

wherein step c is responsive to step a. 
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15. The computer program of Claim 14, wherein step a is calculated as follows: 



llrR 1 . = < 
ml 



- mini jfl, otherwise. 



16. . The computer program of Claim 14, wherein step c is calculated as follows: 

5 

17. The computer program of Claim 14, wherein step c is initialized with soft 
channel information. 

18. The computer program of Claim 14, further comprising the step of providing 
indices of each of the parity check equations. 

10 19. The computer program of Claim 14, further comprising the step of providing an 
index of the smallest LlrQim. 

20. The computer program of Claim 14, further comprising the step of scaling LLrRmi, 
calculated in step a. 

21. The computer program of Claim 13, further comprising the steps of: 
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g converting | LLrAPPk | information into hard information b c ^ 



h calculating an equation syndrome for each parity check equation, Si. 



22. The computer program of Claim 21, further comprising the steps of: 



1 



summing each of equation syndrome in step h, 



5 



J 



outputting the hard information if the sum in step i is equal to zero; and 



k repeating steps a-c if the sum in step i is not equal to zero. 



23. The computer program of Claim 9, further comprising the steps of: 



1 determining if i is less than a first predetermined value; 

m repeating steps a-c, if in step 1, i is less than the first predetermined 

10 value; 



o outputting the hard information if ^ is at least the first predetermined and 
the sum in step i is less than the second predetermined value. 

15 24. The computer program of Claim 23, further comprising the steps of: 

p for each data k, determining if | LLrAPPk | is less than a third 
predetermined value; 

q for each data k> outputting hard information b c> k if | LLrAPPk | is at least 
the third predetermined value; and 

20 r for each data k, outputting soft channel information b s ,k if | LLrAPPk | is 

less than the third predetermined value. 

25. The computer program of Claim 23, further comprising the steps of: 



n 



determining if the sum in step i is less than a second 



predetermined 



value; and 
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s for each data k, determining if a corresponding parity check equation is 
violated; 

t for each data k, outputting hard information b c ,k, if the corresponding 
parity check equation is not violated; 

5 u for each data k } determining if | llrAPPk | is less than a third 

predetermined value; 

v for each data k, outputting hard information b c ,k if | llrAPPk | is at least 
the third predetermined value; and 

w for each data k, outputting soft channel information b S) k if (llrAPPk) is 
10 less than the third predetermined value. 

26. A decoder for decoding low-density parity-check codes comprising: 

first calculating means for calculating UrRmi, for each parity check equation, at 
iteration 

decision aided equalizing means for equalizing LlrRmi, at iteration i, in response 
15 to said first calculating means; and 

second calculating means for calculating LlrQim, for each parity check equation, 
at iteration i in response to said decision aided equalizing means; 

wherein UrQim represents information from bit node 1 to equation node m, one 
for each connection, 

20 wherein UrRmi represents information from equation node m to bit node 1, one 

for each connection, and 

wherein said first calculating means is responsive to said second calculating 

means. 
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27. The decoder of Claim 26, further comprising: 

memory means for storing for each parity check equation, 

the smallest LlrQim calculated by said second calculating means, at 

iteration i; 

the second smallest llrQim calculated by said second calculating means, at 
iteration i\ and 

an overall. sign of LLrQim, calculated by said second calculating means, at 

iteration i; 

wherein said first calculating means is responsive to said memory means, and 
wherein said second calculating means is responsive to said first calculating 

means. 

28. The decoder of Claim 27, wherein said first calculating means calculates as 
follows: 



~ s m 1 • ^{lird \m) • min 1 ' " 1 m , if / * 1^ 1 
1 • sign^tlrQ 1 '/mj ■ mini* ~ * m . otherwise. 



29. The decoder of Claim 27, wherein said second calculating means calculates as 
follows: 
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30. The decoder of Claim 27, wherein said second calculating means is initialized 
with soft channel information. 

31. The decoder of Claim 27, further comprising an address generating means for 
providing indices of each of the parity check equations. 

32. The decoder of Claim 27, wherein said second calculating means provides an 
index of the smallest LlrQim. 

33. The decoder of Claim 27, further comprising multiplying means for scaling 
UrRmi, calculated by said first calculating means. 

34. The decoder of Claim 1, further comprising 

slicing means for converting | llrAPPk | information into hard information b c ,k ; 

and 

equation vector means for calculating an equation syndrome for each parity 
check equation Sj. 

35. The decoder of Claim 34, further comprising: 

summing means for summing each equation syndrome calculated by said 
equation vector means; 

wherein the hard information is output if the sum by said summing means is 
equal to zero, and 
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wherein the calculations by said first and second calculating means are 
repeated if the sum summed by said summing means is not equal to zero. 

36. The decoder of Claim 35, further comprising: 

first threshold means for determining if i is less than a first predetermined 

5 value, 

wherein the calculations by said first and second calculating means are 
repeated if i is less than the first predetermined value as determined by said first 
threshold means; and 

second threshold means for determining if the sum by said summing means is 
10 less than a second predetermined value; and 

wherein the hard information is output if i is at least the first predetermined 
and the sum by said summing means is less than the second predetermined value. 

37. The decoder of Claim 36, further comprising: 

third threshold means for determining for each data k if j LlrAPPk | is less than 
15 a third predetermined value; 

wherein for each data k, outputting hard information b c ,k if | LlrAPPk | is at 
least the third predetermined value as determined by said third threshold means; and 

wherein for each data k, outputting soft channel information b s ,k if | LlrAPPk | is 
less than the third predetermined value as determined by said third threshold means. 

20 38. The decoder of Claim 36, further comprising: 

judging means for determining for each data k if a corresponding parity check 
equation is violated; 

wherein for each data k, outputting hard information b c ,k, if the corresponding 
parity check equation is not violated as determined by said judging means; 
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third threshold means for determining for each data k if | llrAPPk | is less than 
a third predetermined value; 

wherein for each data k, outputting hard information b c ,k if | llrAPPk | is at 
least the third predetermined value as determined by said third threshold means; and 

wherein for each data k, outputting soft channel information b s ^if | llrAPPk | is 
less than the third predetermined value as determined by said third threshold means. 

39. A decoder for decoding low-density parity-check codes comprising: 

a first calculator to calculate UrRmi, for each parity check equation, at iteration 

£-1; 

a decision aided equalizer to equalize UrRmi, at iteration i, in response to said 
first calculator; and 

a second calculator "to calculate UrQim, for each parity check equation, at 
iteration i in response to said decision aided equalizer, 

wherein UrQim represents information from bit node 1 to equation node m, one 
for each connection, 

v wherein UrRmi represents information from equation node m to bit node 1, one 
for each connection, and 

wherein said first calculator is responsive to said second calculator. 

40. The decoder of Claim 39, further comprising: 

a memory to store for each parity check equation, 

the smallest UrQim calculated by said second calculator, at iteration i; 
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the second smallest llrQim calculated by said second calculator, at 
iteration i\ and 

an overall sign of LlrQim, calculated by said second calculator, at iteration 

i\ 

wherein said first calculator is responsive to said memory, and 
wherein said second calculator is responsive to said first calculator. 

41. The decoder of Claim 40, wherein said first calculator calculates as follows: 



IhR 



i - 
ml 



1 • sign[llrd 1 /„) • min 1 ' " 1 m , if / * / 



j-1 
m 



-i 1 ■ 



(llrQ 1 l im) 



■ mini' * m , otherwise. 



42. The decoder of Claim 40, wherein said second calculator calculates as follows: 

UrQj-UrP}+ £ 

43. The decoder of Claim 40, wherein said second calculator is initialized with soft 
channel information. 
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44. The decoder of Claim 40, further comprising an address generator for providing 
indices of each of the parity check equations. 

45. The decoder of Claim 40, wherein said second calculator provides an index of the 
smallest LlrQun. 

5 46. The decoder of Claim 40, further comprising a multiplier to scale LLrRmi, 
calculated by said first calculator. 

47. The decoder of Claim 39, further comprising 

a slicer to convert | LlrAPPk | information into hard information b Ct k ; 

an equation vector circuit to calculate an equation syndrome for each parity 
10 check equation, Si. 

48. The decoder of Claim 47, further comprising: 

a summer to sum each equation syndrome calculated by said equation vector 

circuit; 

wherein the hard information is output if the sum summed by said summer i is 
15 equal to zero; and 

wherein the calculations by said first and. second calculator are repeated if the 
sum summed by said summer is not equal to zero. 

49. The decoder of Claim 48, further comprising: 

a first threshold detector to determine if i is less than a first predetermined 

20 value; 

.wherein the calculations by said first and second calculator are repeated if i is 
less than the first predetermined value as determined by said first threshold detector; 
and 
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a second threshold detector to determine if the sum by said summer is less 
than a second predetermined value; and 

wherein the hard information is output if i is at least the first predetermined 
and the sum by said summing means is less than the second predetermined value. 

5 50. The decoder of Claim 49, further comprising: 

a third threshold detector to determine for each data k if | llrAPPk | is less than 
- a third predetermined value; 

wherein for each data k } outputting hard information b c ,k if | llrAPPk | is at 
least the third predetermined value as determined by said third threshold detector; 
10 and 

wherein for each data k, outputting soft channel information b s ,k if | llrAPPk | is 
less than the third predetermined value as determined by said third threshold 
detector. 

51. The decoder of Claim 49, further comprising: 

15 a judging circuit to determine for each data k if a corresponding parity check 

equation is violated; 

wherein for each data k, outputting hard information b c ,k, if the corresponding 
parity check equation is not violated as determined by said judging circuit; 

a third threshold detector to determine for each data k if | llrAPPk | is less than 
20 a third predetermined value; 

wherein for each data k, outputting hard information b c ,k if | llrAPPk | is at 
least the third predetermined value as determined by said third threshold detector; 
and 
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wherein for each data k, outputting soft channel information b s ,k if | LlrAPPk | is 
less than the third predetermined value as determined by said third threshold 
detector. 

52. The method of Claim 1, further comprising the steps of: 

g converting | LlrAPPk | information into hard information b c ,k; 

h calculating an equation syndrome for each parity check equation, Si. 
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ABSTRACT 

A. decoder for decoding low-density parity-check codes is provided comprising a 
first calculator to calculate LlrRmi, for each parity check equation, at iteration a 
decision aided equalizer to equalize UrRmi, at iteration i, in response to the first 
5 calculator; and a second calculator to calculate llrQim, for each parity check equation, 
at iteration i in response to the decision aided equalizer. llrQim represents information 
from bit node 1 to equation node m, one for each connection, and UrRmi represents 
information from equation node m to bit node 1, one for each connection. The first 
calculator is responsive to the second calculator. The decoder further comprises a 
10 memory to store for each parity check equation, the smallest llrQim calculated by the 
second calculator, at iteration i; the second smallest LlrQim calculated by the second 
calculator, at iteration i; and an overall sign of llrQim, calculated by the second 
calculator, at iteration i. The first calculator is responsive to the memory, and the 
second calculator is responsive to the first calculator. 
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